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Understanding and controlling the way excitons propagate in solids is a key for tailoring materials 
with improved optoelectronic properties. A fundamental step in this direction is the determination of 
the exciton energy-momentum dispersion. Here, thanks to the solution of the parameter-free Bethe- 
Salpeter equation (BSE), we draw and explain the exciton energy-momentum map of hexagonal 
boron nitride (h-BN) in the hrst three Brillouin zones. We show that h-BN displays strong excitonic 
effects not only in the optical spectra at vanishing momentum q, as previously reported, but also at 
large q. We validate our theoretical predictions by assessing the calculated exciton map by means 
of an inelastic x-ray scattering (IXS) experiment. Moreover, we solve the discrepancies between 
previous experimental data and calculations, proving then that the BSE is highly accurate through 
the whole momentum range. Therefore, these results put forward the combination BSE and IXS as 
the tool of choice for addressing the exciton dynamics in complex materials. 


Introduction - The response of materials to electro¬ 
magnetic fields is determined by electronic excitations 
that are strongly influenced by electron-hole (e-h) inter¬ 
actions. In particular, the e-h attraction leads to the 
formation of excitons, which are a fundamental aspect 
in the functionality of many optoelectronic devices, as 
excitons can propagate in materials carrying excitation 
energy that can be transformed and exploited by different 
means. Excitons can be identified in electronic spectra 
as sharp peaks within the band gap of insulators mm or, 
beyond the band gap, as spectral intensity enhancement 
towards lower energies with respect to a non-interacting 
theoretical picture [3]. A reliable description and analy¬ 
sis of those two-particle correlation effects is therefore the 
key to understand materials’ dielectric properties, guide 
the realisation of new experiments and foster the devel¬ 
opment of new technological applications. 

Nowadays, the ab initio solution of Bethe-Salpeter 
equation (BSE) |4HZ] represents the state-of-the-art 
method to obtain spectra in very good agreement with 
experiments in a large variety of materials [8]. These 
theoretical achievements have mainly focused on optical 
absorption and electron energy-loss spectroscopy (EELS) 
spectra for vanishing momentum transfer q ^ 0. How¬ 
ever, judging whether a theoretical approach captures 
fully the physics of the electron dynamics requires strin¬ 
gent tests that are offered by the measurement of the 
full dynamic range of the relevant variables (momentum 
and energy). This assessment is now possible thanks 
to the spectacular progress of inelastic x-ray scattering 
(IXS) experiments, both in the resonant (RIXS) and non¬ 
resonant (NRIXS) conditions. They allow one to probe 
electronic excitations at finite momenta (q) [9] with a re¬ 


solving power that has improved by orders of magnitude 
in the last two decades HDHH]. The challenge for theory 
is hence the first-principles description of the full electron 
dynamics m, which is well beyond the sole simulation 
of optical absorption spectra. The investigation of the 
energy-momentum dispersion of elementary excitations 
(excitons, plasmons, etc.) provides fundamental informa¬ 
tion on the way they propagate in materials. Moreover, 
spectroscopic features measured at larger q allow for the 
study of excitations that in real space occur on shorter 
interatomic scales and/or are not visible in optics because 
they are dipole forbidden (see e.g. Refs. pMHIfi] ). 

Recently, some of us have shown that these ambitious 
goals are within reach, thanks to the extension of the BSE 
to describe e-h excitations carrying a finite momentum 
q m- In the case of a prototypical wide-gap insulator 
such as lithium fluoride, the dynamic structure factor 
5'(q, cj) was obtained in excellent agreement with accu¬ 
rate NRIXS data [18]. However, in contrast to LiF, in 
hexagonal boron nitride (h-BN), a layered material that 
is the insulating counterpart of graphite, a recent com¬ 
parison m between BSE calculations and NRIXS data 
revealed a mismatch at high momentum transfers. The 
experiment displays a sharp peak at ^ 7 eV whose ori¬ 
gin has still to be understood, as pointed out in Ref. [191 
Moreover, despite the accurate BSE results in LiF HI], a 
recent study based on a simplified exciton kinetic kernel 
model [20] raised doubts about the capability of the ab 
initio BSE in general to address the issue of the exciton 
band structure. Therefore new questions arise: Is LiF 
only a fortunate case for BSE? For spectra at finite q 
should we hence go beyond the standard BSE implemen¬ 
tation that has been successfully applied for optics? Or, 
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on the contrary, is the sharp peak measured in h-BN an 
artifact of the NRIXS experiment? Finally, is the BSE 
the appropriate method in general to study the exciton 
band structure of complex materials? 

In the present work, we demonstrate that it is possi¬ 
ble to reconcile theory and experiment if the microscopic 
details of the screened Coulomb e-h interaction are ex¬ 
plicitly taken into account in the calculations. We show 
that in h-BN excitonic effects not only produce a redshift 
of the peaks, as observed in Ref. |T9l but also lead to im¬ 
portant spectral shape redistributions. We reproduce all 
the details of the NRIXS spectra measured in Ref. [191 
including the 7 eV peak at large q, and we consistently 
explain the appearance of excitonic features at various q 
as a multiple manifestation of the same large joint den¬ 
sity of states (JDOS). Moreover, by calculating the entire 
map of 5'(q, cj) in the first 3 Brillouin zones of h-BN, we 
identify new “hot spots” in dynamics, where the proba¬ 
bility for creating an exciton is the largest. These predic¬ 
tions are fully confirmed by the new NRIXS experiment 
that we have performed at the beamline ID20 of the Euro¬ 
pean Synchrotron Radiation Eacility (ESRE) in Grenoble 
(Erance). Therefore these new results demonstrate that 
the combination of first-principles BSE calculations and 
accurate NRIXS experiments is a very powerful tool to 
explore and understand the exciton dynamics. 

Theory - The spectrum obtained in an NRIXS exper¬ 
iment is proportional to the dynamic structure factor: 


5'(q,w) 


(47r2n) 




( 1 ) 


where n is the average electron density. The inverse of 
the macroscopic dielectric function e'^ can be expressed 

as [HI[TT]: 


— 1 / \ -1 
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where in the Tamm-Dancoff approximation (TDA) [21] 
the sum is over valence-conduction {v-c) transitions t, 
the oscillator strengths are p^(q) = (0vk-qr 
k and qr are in the first Brillouin zone, and q = 
qr + G is the measured momentum transfer with a 
reciprocal-lattice vector G. The BSE can be cast into 
an effective two-particle Schrodinger equation H E]: 
iJexc(q)-4A(q) = -BA(q)AA(q), where ^A(q) and £’A(q) 
are the exciton eigenvectors and eigenvalues, respec¬ 
tively. The exciton hamiltonian i^exc- (^l^excl^O ~ 
E {t\vc — W\t') contains the quasiparticle e-h tran¬ 
sition energies Et calculated in the GW approximation 
[22] , and the matrix elements in the transition basis of the 
bare Coulomb interaction Vc and the statically screened 
Coulomb interaction W = €~^Vc, which describes the e-h 
attraction (i.e. gives rise to excitons) and is here obtained 


in the random phase approximation (RPA). We have cal¬ 
culated the ground state of h-BN within the local density 
approximation (EDA) [23| of density functional theory 
[24] . using norm-conserving Troulliers-Martin pseudopo¬ 
tentials [25] in a plane wave approach [26] with an energy 
cutoff of 30 Hartree. In order to approximate GW quasi¬ 
particle energies, following m we corrected the LDA 
band structure with a scissor operator of 1.98 eV for the 
band gap and stretching the valence bands by 5% . The 
BSE spectra at finite q are obtained from the EXC code 
[28] using 25 bands and a 18 x 18 x 4 k-point grid for 
calculations with q parallel to the hexagonal plane and 
with al2xl2x8 k-point grid for calculations with q 
perpendicular to it. 

Experiment - In the NRIXS experiment the incident 
x-ray beam from three undulators was monochromated 
to an energy of 7.5 keV by a combination of a Si(lll) 
double crystal and a Si(311) channel-cut. The beam was 
focused to a spot of « 10 /im x 20 /im (V x H). The 
spectrometer used a diced Si(533) analyzer crystal in the 
Johann geometry with a Rowland circle diameter of 2 
m. The active diameter of the analyzer crystal was 80 
mm, yielding a momentum-transfer resolution of ^ 0.15 
A“^. The energy resolution was 200 meV (EWHM). The 
scattering plane was vertical, i.e., perpendicular to the 
plane of linear polarization of the incident and scattered 
radiation. The detector was based on a Timepix chip, 
with a pixel size of 55 /im and enabling us to use the dis¬ 
persion compensation algorithm [29l |30] . The analyzer 
Bragg angle was fixed at 87° and the energy-loss spectra 
were measured by scanning the incident-photon energy. 
The h-BN single crystal was a colorless and transparent 
platelet with a thickness of 0.5 mm and lateral size of 2.1 
mm, which was produced by a high- pressure and high- 
temperature (HPHT) method using a barium-related sol¬ 
vent system as reported in m- The samples were aligned 
using x-ray diffraction on the beamline. 

Results - In layered h-BN the electronic states can be 
classified according to the even (a) and odd (tt) parity 
with respect to the single BN sheet [32]. Optical and 
EELS spectra [33ll36] in the low-energy range, where ex¬ 
citonic effects are more relevant, are determined by tt—tt* 
transitions that are mainly visible for in-plane q —> 0. 
The weak screening of the Coulomb e-h interaction, in 
concomitance with a JDOS peak due to vertical transi¬ 
tions between parallel bands m in the HKML plane of 
the Brillouin zone [32l [38] , in the q ^ 0 spectra gives rise 
to prominent exciton peaks with large binding energies 

EaEiiio]. 

Fig. 0 shows the dynamic structure factor 5'(q, cj) in 
the energy range of the tt — tt* transitions at three finite 
momentum transfers: q = A, K, and M, which are lo¬ 
cated at the boundaries of the first Brillouin zone m- 
As noticed in m, the e-h attraction induces a redshift 
of the peaks, which is larger for the in-plane q direc¬ 
tions at K and M (redshift of 1.5 eV) than for q at A 
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FIG. 1. Dynamic structure factor at three Brillouin zone 
boundaries q = A, K, and M calculated in the GW-RPA and 
from the solution of the BSE compared to the NRIXS data 
from Ref. [19]. 


(redshift of 0.4 eV), i.e. along the direction perpendic¬ 
ular to the BN layers, implying an anisotropic effect of 
the e-h interaction m- This is evident from the com¬ 
parison of the BSE spectra with the GW-RPA results, 
obtained starting from the GW band structure and ne¬ 
glecting the e-h attraction W in the exciton hamiltonian 
ii^exc- At variance with m, where the shift was inferred 
from the adjustment of the calculations to the experimen¬ 
tal spectra, the anisotropic redshift of the spectra is here 
the direct outcome of the BSE calculations that result 
in very good agreement with NRIXS. Our simulations of 
the experimental spectra, being free from any adjustable 
parameters, allow us to additionally observe that exci- 
tonic effects also induce an important redistribution of 
the spectral weight towards lower energies with respect 
to the non-interacting e-h picture. 

At q = M the two main peaks are located at ^ 8 and 
^12 eV. At the K point, the first peak at ^ 6-7 eV 
appears as a shoulder of the second one. In both cases 
the 12 eV peak originates from non-vertical (i.e. k ^ 
k' = k + q) TT — TT* transitions that disperse isotropically 
as a function of in-plane q, shifting the energy of the 
peak from 9 eV at q ^ 0 to 12 eV in both EK and EM 
directions. The first peak, instead, is due to a peculiar 
property of the hexagonal Brillouin zone, see Fig. |^a). 
For q =M it derives from vertical e-h transitions between 
k points belonging to the ML line in the band structure 
[37] (analogous, for q=K, is the HK line). These are the 
same vertical tt —tt* transitions with large JDOS that are 
at the origin of the tightly bound exciton in the spectra 
at q ^ 0 [27[ [39| . The same strong excitonic effects are 
hence appearing also at q =M and K as a spectral shape 
redistribution that similarly strongly enhances the first 
peak in the spectrum. They create a new shoulder at 
q =K (entirely absent in GW-RPA) and at q =M they 
make the first peak become more intense than the second 
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FIG. 2. Dynamic structure factor for q along FK direction 
calculated in the GW-RPA and from the solution of the BSE 
compared to the NRIXS data from Ref. [19]. The q = 2.5 
is the M' point (see text). 


one, in contrast to the GW-RPA results. 

Having clarified the diverse role of excitonic effects at 
various q, we are ready to discuss what happens along 
the FK direction, where for q > 2.0 a peak appears 
in the NRIXS spectra at ^ 7 eV that was unexplained 
in Ref. [19] . The spectra that we have calculated for dif¬ 
ferent momentum transfers along FK demonstrate that 
the BSE is actually able to reproduce the 7 eV peak that 
is measured in the NRIXS experiment US] (see Fig. |^. 
Moreover they show that the peak has indeed an exci¬ 
tonic character as it is completely absent in the GW-RPA 
calculations. In order to better understand the origin of 
this excitonic peak one has to note that it has the largest 
intensity at q = 2.5 A“^. This momentum transfer along 
the FK direction is in fact another M point, which is lo¬ 
cated at the boundary of the third Brillouin zone and 
which we call M'. We can also immediately recognise 
why excitonic effects are again particularly strong at this 
momentum transfer M'. The exciton eigenvalues Ex{q) 
and eigenvectors Ax (q) at M' must be the same as at the 
M point in the first Brillouin zone since the two points 
just differ by a reciprocal-lattice vector G. However the 
oscillator strengths p'^{q) are generally different at dif¬ 
ferent q and this explains why the spectra at M and M' 
do not entirely overlap. In particular, at q=M' there is 
a prominent peak at 7 eV emerging from a featureless 
plateau at higher energy, while at q=M a double-peak 
structure was observed in the same energy range. 

We have therefore found that the solution of the BSE, 
when the screening of the Coulomb interaction is explic¬ 
itly calculated at the RPA level rather than obtained 
from model dielectric functions dnumm, is able to re- 
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FIG. 3. (a) Color map of S'(q,ct;) calculated from the BSE for fixed cj = 7 eV and for all q spanning the first 3 Brillouin zones 
with qz = 0. To enhance its visibility, AS'(q, c^;) in the hrst BZ has been magnihed by a factor 3. The high-symmetry points 
relevant for the discussion in the main text are explicitly indicated (in parenthesis those belonging to the q^ = 0.5 plane). 
Adding q = TM (solid violet line) to k =M (L) one obtains another k' = k + q =M (L) (see dashed violet line). Equivalently 
adding q = TK (solid red line) to k =K (H) one obtains again k' =K (H) (see red dashed line), (b) Comparison between 
calculated and measured *S'(q,Ct;) for q along the TM direction. 


produce and explain the NRIXS data. This is confirmed 
by the direct comparison between the calculated and ex¬ 
perimental spectra m- The agreement is excellent in 
all considered cases. Moreover, having understood that 
strong excitonic effects for the peak at 7 eV should be 
expected whenever tt — tt* transitions give rise to an in¬ 
tense JDOS, we can forecast other “hot spots” with large 
intensity in the dynamic structure factor, beyond what 
has been experimentally detected in m- In Fig. [^a) we 
thus plot a color map of 5'(q, uj) for a fixed uj = 7 eV and 
for all q spanning the first 3 q;^ = 0 Brillouin zones. In 
this manner we can easily identify two other remarkable 
points along the FM direction, located between q=M and 
the vertex of the hexagon q=F', which differs from the F 
point by a reciprocal-lattice vector G. 

In order to confirm these predictions and assess the 
whole theoretical map of 5'(q, cj) that has been obtained 
from the BSE calculations, we compare the theoretical 
results to those obtained from our NRIXS experiment 
in Fig. I^b). While the GW-RPA results are totally 
different from the experiment, the excellent agreement 
between the BSE results and the NRIXS data proves the 
predictive power of the BSE. In the spectra the largest 
intensity is found at q = F', displaying two prominent 
peaks deriving again from the large JDOS of vertical 
TT — TT* transitions. The first peak, visible only close to 
q = F', matches the main exciton peak at the onset of 
the spectrum at q ^ 0. The second peak at q = F', 
which is not noticeable at q ^ 0, instead evolves contin¬ 
uously from the first peak located at ^ 8 eV at q=M. 
These results show that inspecting spectra at large q, 
thanks to the variation of the oscillator strengths with q, 
can reveal formation of excitons that are hidden in corre¬ 
sponding optical spectra. Moreover, they evidence how 


the combination of NRIXS and BSE provides the means 
for their detection in the whole energy-momentum range, 
allowing one to obtain the full exciton band structure. 

Conclusions - In summary, by solving the discrepan¬ 
cies between previous experiments and calculations, we 
have successfully established the energy-momentum map 
covering the first 3 Brillouin zones of the dynamic struc¬ 
ture factor 5'(q, cj) of h-BN, a prototypical layered in¬ 
sulator. We have shown that excitonic effects in h-BN 
are strong also at large momentum q. They are essen¬ 
tial to interpret and understand the spectra (including 
previous unexplained features) and cannot be neglected 
remaining at a level of theory corresponding to the RPA 
m- We have explicitly proved that such a theoretical 
map based on BSE is crucial in order to guide the exper¬ 
imental exploration of the electronic dynamics, as new 
IXS measurements have fully validated the theoretical 
predictions. Therefore these case-study results promote 
the ab initio solution of the BSE, which has been already 
successfully applied in optics, as an accurate and predic¬ 
tive method also to investigate the charge dynamics and 
the exciton band structure for a wide range of materials 
and technological applications. 
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